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ABSTRACT 

Context. Dust attenuation curves in external galaxies are useful to study their dust properties as well as to interpret their intrinsic 
spectral energy distributions. These functions are not very well known in the UV range either at low or high redshift. In particular 
the presence or absence of a UV bump at 2175 A remains an open issue which has consequences on the interpretation of broad band 
colours of galaxies involving the UV range. 

Aims. We want to study the dust attenuation curve in the UV range at z > 1 where the UV is redshifted into the visible, and with 
Herschel data to constrain dust emission and a global dust attenuation. In particular we will search for the presence of a UV bump 
and the related implications on dust attenuation determinations. 

Methods. We use deep photometric data of the Chandra Deep Field South obtained with intermediate and broad band filters by the 
MUSYC project to sample the UV rest-frame of galaxies with 1 < z < 2. Herschel/FACS and Spitzer/MIFS data are used to measure 
the dust emission. 30 galaxies were selected with high S/N in all bands. Their SEDs from the UV to the far-IR are fitted using the 
CIGALE code and the characteristics of the dust attenuation curve are obtained as Bayesian outputs of the SED fitting process. 
Results. The mean dust attenuation curve we derive exhibits a significant UV bump at 2175 A whose amplitude corresponds to 
35% (76 %) that of the Milky Way (Large Magellanic Cloud: LMC2 supershell) extinction curve. An analytical expression of the 
average attenuation curve (A(A)/Ay) is given, it is found slightly steeper than the Calzetti et al. one, although at a Icr level. Our galaxy 
sample is used to study the derivation of the slopes of the UV continuum from broad band colours, including the rest-frame GALEX 
FUV - NUV colour. Systematic errors induced by the presence of the bump are quantified. We compare dust attenuation factors 
measured with CIGALE to the slope of the UV continuum and find that there is a large scatter around the relation valid for local 
starbursts (~ 0.7 mag). The uncertainties on the determination of the UV slope lead to an extra systematic error of the order of 0.3 to 
0.7 mag on dust attenuation when a filter overlaps the UV bump. 
Conclusions. 

Key words. Galaxies: high-redshift-Galaxies: ISM-Ultraviolet : galaxies-Infrared: galaxies-ISM: dust, extinction 






1. Introduction 

Stellar light in galaxies is absorbed and scattered by the interstel- 
lar medium, because of the presence of dust grains. On a galaxy 
scale, the process is usually quantified by the introduction of an 
average attenuation function which reflects the combined eflTects 
of absorption and scattering in a complex geometry. 
The most important feature in the extinction curves of the Milky 
Way (MW) and Large Magellanic Cloud (LMC) is the so called 
UV bump, a strong absorption feature at 2175 A which is not 
observed in the Small Magellanic Cloud (SMC) bar, although a 
line of sight through the SMC wing exhibits an extinction curve 



with a strong UV bump (Li & Draine 2002 Cartledge et al. 



|2005l). The e xact origin of this feature is still under discussion 
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(e.g. |Drailie| |1989[ |Xiang etaT, "2011" and references therein) 
When dealing with average attenuation curves, the search for 
a bump is complicated by the eflTects of scattering or geometry 
dlnoue et al.| 120061 |NoU et aL|[2007l|Panuzzo et al.||2007 ) which 
may aflPect the amplitude of the resulting bump in an average 
attenuation curve even when assuming a local Milky Way like 
extinction curve. 

The best way to study UV bumps in external galaxies is to per- 
form UV rest-frame spectroscopy: Calzetti et al. ( 1994| ) deduced 
the net attenuation curves of local starburst galaxies by compar- 
ing the UV spectra of dusty and dust- free spe ctra of local star- 
bursts: they did not find any bump. Conversely, |Noll etaL (2005 
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2009a| observed a significant 2175 A absorption feature in the 
spectra of star forming galaxies at z ~ 2 consistent with an ex- 
tinction curve intermediate between the SMC and LMC ones. 
The 2175 A feature has also been detected in the spectral energy 
distribution of some distant galaxies hosting gamma-ray bursts 
( [Liang & Lil|2009[[2010| ), the exti nction curves are found to dis- 



3200 



play a large diversity of shapes (L iang & Li| |2QQ9| |2QlQ[|Zafar 
etal.||2Qll and references therein). We also refer to Xiang et al. 



] 2011| ) foTan extensive review on the detection of the UV bump 
in Galactic and extragalactic interstellar medium. 
Broad band colours are also used to search for evidence of a UV 
bump. From an analysis of the observed B-R colours of galaxies 
at z ~ 1 Conroy (2010) disfavours the pr esence of UV bumps as 
strong as that in the Milky Way whereas |Burgarella et al. ( 2005| ) 
analysed the broad bands SEDs of nearby galaxies observed by 
GALEX and IRAS and conclude to the presence of a bump in 
the attenuation curve of these galaxies. In all cases, the use of 
broad band colours which overlap the bump makes difficult to 
disentangle the eff'ects of a dust extinction curve with or with- 
out a bump, of dust/stars geometry and of various star forma- 
tion histories ([Inoue et al.> 2006 ; Panuz zo et al.[ [2007 [ |Johnson| 



et al. 2007a^bi Buat et al. 



2011| ) leading to contradicting con- 



clusions on the presence of a UV bump in nearby galaxies ob 
ser ved by G ALEX and SDSS ( [Conroy et aLJ [20Tq1 [Johnson et 
|aL| [2007a|b| ). 

The presence and characteristics of the UV bump, and far UV 
rise in the extinction and attenuation curves can provide infor- 
mation on the interstellar medium properties ( [Noll et aL|[2009a| . 
Beyond the studies of the interstellar medium, an reliable atten- 
uation curve in the UV range is very useful to derive star forma- 
tion rates in galaxies since the UV emission is commonly used 
as a star formation tracer. For example, [Wij esinghe et al. [ pO 1 ) 
compared star formation rates deduced from luminosities mea- 
sured in the GALEX bands and the Ha and [Oil] emission lines 
and concluded that they must remove the 2175 A feature from 
the attenuation curve to obtain the best agreemen t between the 
diff'erent SFR estimators. The attenuation curve of lCalzetti et al.l 
( [2000] ) is commonly used, especially in high-z studies and does 
not include a bump: if such a bump is present, omitting it leads 
to incorrect dust attenuation corrections (Ilbert et al. , 2009 ). 
The global amount of dust attenuation is robustly estimated by 
comparing dust and stellar emission, through the Lir/Luv ra- 
tio. When IR data are not available the shape of the UV con- 
tinuum (< 3000 A) has b een proposed as a p roxy to measure 
the amount of attenuation. [Calzetti et al.[ ( p^94| ) showed that the 
UV continuum of local starburst galaxies can be fit by a power 
law (/i(erg cm-^s-^A"^) oc A^) for A > 1200 A. They calculated 
13 by defining 10 windows in the lUE spectra of th eir galaxies 
avoiding the 2175 A feature. Meurer et al. (1999]) found a re- 
lation between jS and the amount of dust attenuation measured 
by Lir/Luv- This local starburst relation is widely used to es- 
timate dust attenuation in high redshift galaxies using a broad 
band colour in the UV rest-frame of the galaxies as a proxy of 



p (e.g [Burgarella et al.| [20071 [Elbaz et al.[ [2007[ [Daddi et al 



2007 [ [Reddy et al.[ [2008[ ). If a bump is present in the dust at- 
tenuation curve, one must be cautious by choosing broad band 
filters which do not overlap the bump in order to estimate reliable 
slopes outside the wavelength range affected by the bump, as was 
done to derive the [Meurer et al.[ ( p^ 9 ) relation. This has been il- 
lustrated by studies of the local universe based on GALEX FUV 
and NUV data: the NUV band largely overlaps the UV bump and 
the FUV-NUV colour is found to be sensitive to the presence or 
absence of this feature as well as to the recent star formation his- 
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Fig. 1: Rest- frame ultraviolet wavelengths sampled with our 
dataset as a function of redshift. The solid curves correspond to 
broad band filters (from bottom left corner to top right corner: 
U,U38,B,V,R,I and z) and the dashed curves to intermediate- 
band filters (from bottom left corner to top right corner: lA- 
427,445,484,505,527,550,574,598,624,651 (coincident with the 
red soHd curve for the R-band), 679 and 738. The FWHM of 
the intermediate filters varies from 20 to 34 nm. The horizontal 
dotted line corresponds to 2175 A, the central wavelength of the 
expected bump. 



tory ( [Panuzzo etaL) [2007 [[Inoue et al.[[2006| ). 
In this paper we take advantage of the availability of high quality 
deep optical data obtained through intermediate-band filters in 
the Chandra Deep Field South (CDFS) to sample the rest-frame 
UV of galaxies at high redshift. These data are combined with 
the Herschel observations of the CDFS field as part of the open 
time key program GOODS-Herschel project (RI. D. Elbaz). The 
galaxies detected in the far-IR are expected to experience a sub- 
stantial extinction, which helps to detect the imprint of the bump 
in their SED. The combination of UV and IR data will allow us 
to compare the properties of the UV continuum to the dust at- 
tenuation at work in these objects. In a similar vein, Ilbert et al. 
(2009) used broad and intermediate band filters to derive pho- 
tometric redshifts in the COSMOS field and found that a UV 
bump is necessary to explain the UV rest-frame fluxes of some 
starburst galaxies but they did not discuss the amplitude of the 
bump and dust attenuation in these galaxies. 
The construction of the sample is explained in section 2. We 
want to use only sources detected with high signal to noise ratio 
(SNR) in a large number of filters which sample well the rest- 
frame UV wavelength region aff'ected by the bump. This cri- 
terium leads to a redshift selection between 0.95 and 2.2. The 
dust attenuation curve as well as the amount of dust obscuration 
is constrained by fitting the whole SEDs of the galaxies with 
the CIGALE code, which is described in section 3. Section 4 is 
devoted to the results about the presence of the UV bump and 
section 5 to the consequences on the determination of the slope 
of the UV continuum and dust attenuation in galaxies. The con- 
clusions are presented in section 6. 

2. Selection of the sample 

As part of the GOODS-Herschel key program ([Elbaz et al.| 



2011), the Herschel Space Observatory (Pilbratt et al. 2010] 



surveyed a small area of the Great Observatories Origins Deep 
Survey Southern field (GOODS-S): 10' x 10' centered on the 
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Table 1 : Galaxy sample 



N 


RA (2000) 
deg 


DEC (2000) 
deg 


z 


log(LiR) 


Auv 
mag 


£b 


errE, 


r 

A 


err^ 

A 


S 


errs 




1 


03:32:36.19 


-27:42:58.2 


0.967 


11.45 


1.14 


0.98 


0.67 


339.32 


107.00 


0.09 


0.03 


3.16 


2 


03:32:20.72 


-27:44:53.7 


0.969 


11.50 


2.80 


1.74 


0.31 


361.66 


102.58 


-0.30 


0.00 


2.69 


3 


03:32:16.45 


-27:44:48.7 


0.976 


11.49 


2.98 


1.15 


0.61 


348.02 


107.33 


0.03 


0.07 


2.63 


4 


03:32:31.80 


-27:49:58.4 


0.98^ 


11.49 


2.89 


1.18 


0.54 


360.33 


107.93 


-0.16 


0.09 


1.43 


5 


03:32:42.28 


-27:47:46.1 


0.998 


11.43 


1.98 


1.68 


0.34 


376.98 


101.85 


-0.20 


0.12 


2.31 


6 


03:32:37.39 


-27:44:07.0 


1.018 


11.28 


3.90 


1.51 


0.46 


350.17 


106.16 


-0.24 


0.06 


1.76 


7 


03:32:06.43 


-27:47:28.7 


1.021 


11.69 


3.75 


1.25 


0.54 


370.61 


106.19 


-0.29 


0.03 


2.30 


8^ 


03:32:15.78 


-27:46:30.2 


1.023 


11.35 


1.68 


1.05 


0.65 


347.67 


107.83 


-0.21 


0.05 


3.72 


9 


03:32:36.02 


-27:44:23.8 


1.038 


11.34 


1.98 


1.55 


0.45 


351.24 


106.58 


-0.02 


0.05 


1.65 


10 


03:32:10.58 


-27:47:06.1 


1.045 


11.42 


2.83 


0.95 


0.62 


352.28 


108.69 


-0.10 


0.09 


1.24 


11 


03:32:28.98 


-27:49:08.2 


1.094 


11.59 


2.01 


1.63 


0.43 


375.24 


102.58 


-0.29 


0.02 


5.58 


12 


03:32:39.87 


-27:47:15.0 


1.095 


11.43 


1.93 


1.79 


0.29 


400.42 


92.73 


-0.27 


0.07 


2.25 


13 


03:32:25.86 


-27:50:19.6 


1.095 


11.52 


2.40 


1.42 


0.52 


355.97 


106.19 


-0.16 


0.06 


2.30 


141 


03:32:17.18 


-27:52:20.5 


1.097 


12.09 


5.61 


1.03 


0.39 


368.57 


104.51 


-0.30 


0.02 


4.56 


15 


03:32:34.85 


-27:46:40.5 


1.099 


11.46 


4.19 


1.61 


0.43 


373.69 


102.00 


-0.18 


0.08 


1.08 


16 


03:32:34.63 


-27:53:24.5 


1.107 


11.61 


2.47 


1.04 


0.64 


344.75 


107.87 


-0.03 


0.10 


1.12 


171 


03:32:29.97 


-27:45:29.8 


1.209 


12.07 


1.13 


0.87 


0.67 


328.56 


106.76 


0.08 


0.04 


15.92 


18 


03:32:43.48 


-27:45:56.5 


1.220 


11.45 


2.54 


1.35 


0.55 


336.36 


105.68 


-0.21 


0.06 


2.70 


19 


03:32:47.17 


-27:51:06.1 


1.224 


11.82 


4.66 


1.10 


0.50 


337.70 


107.00 


-0.29 


0.03 


1.22 


20^ 


03:32:35.97 


-27:48:50.4 


1.309 


11.54 


1.24 


0.87 


0.67 


337.88 


107.68 


0.02 


0.07 


2.70 


211 


03:32:29.93 


-27:43:01.0 


1.372 


12.12 


4.29 


1.62 


0.32 


408.06 


89.38 


0.01 


0.06 


1.71 


22^ 


03:32:28.79 


-27:47:55.6 


1.383 


11.69 


4.02 


0.66 


0.47 


331.44 


107.91 


-0.21 


0.07 


3.36 


23 


03:32:34.03 


-27:50:28.8 


1.384 


11.82 


3.32 


1.45 


0.38 


369.56 


102.64 


0.08 


0.04 


2.50 


24 


03:32:41.33 


-27:45:38.2 


1.522 


11.91 


3.69 


1.30 


0.51 


332.05 


107.31 


-0.20 


0.08 


1.53 


251 


03:32:37.76 


-27:52:12.3 


1.609 


12.33 


4.63 


1.03 


0.51 


333.88 


111.60 


-0.04 


0.05 


6.39 


26 


03:32:37.71 


-27:50:00.4 


1.619 


12.43 


5.21 


0.94 


0.36 


352.65 


105.73 


-0.09 


0.07 


2.88 


27 


03:32:17.82 


-27:46:39.7 


1.956 


11.93 


2.85 


0.83 


0.60 


368.08 


107.77 


0.01 


0.08 


1.69 


28 


03:32:41.77 


-27:46:56.4 


1.994 


12.00 


3.32 


1.03 


0.50 


351.00 


107.88 


0.02 


0.09 


2.85 


29 


03:32:40.05 


-27:47:55.4 


1.997 


12.40 


3.79 


1.39 


0.40 


328.89 


103.03 


-0.13 


0.10 


2.30 


30 


03:32:40.75 


-27:49:26.1 


2.130 


12.12 


3.02 


1.33 


0.53 


350.01 


107.06 


-0.08 


0.11 


1.41 



Notes. IR luminosities(log(LiR)), dust attenuation (Auv) and parameters relative to the attenuation curve (the amplitude and width of the UV 
bump,£^^ and y, and the steepness of the curve, 6) and their corresponding errors are obtained from the Bayesian analysis performed with CIGALE. 
^^^ X-ray emitter 
*^2^ No spectroscopic redshift, the photometric redshift is used. 



CDFS were obser ved at 100 and 1 60 //m over 264 hours by the 
PACS instrument ( Poglitsch et al. |2010^ . Source extraction on 
the PACS images was performed at the prior position of Spitzer 
24 yum sources as described in |Elbaz et al.| ( [201 \\ . 
Cardamone et al. (2010, the MUSYC project) compiled a uni- 
form catalog of optical and infrared photometry for sources in 
GOO DS-S, incorporating the G OODS Spitzer IRAC and MIPS 
data ( [Dickinson et al.| 2003| ). In particular, they used deep 
intermediate-band imaging from the Subaru telescope to pro- 
vide photometry with finer wavelength sampling than is possible 
from standard broad band data. This was done in part to enable 
more accurate photometric redshifts, but for our purposes here 
it provides a valuable means of tracing the detailed shape of the 
UV rest-frame spectrum and to look for the presence of a UV 
absorption bump. 

We start with the GOODS-Herschel catalogue of 533 sources 
detected at 100 //m, and identify 325 sources with a SNR (flux 
over error) at 100 yum larger than 5. We further restrict the sam- 
ple to the 235 sources also detected at 160 yum with a SNR 
> 3 and we cross-correlate them with the MUSYC catalogue of 
84402 sources with photometric redshifts. The tolerance radius 
between IRAC and optical coordinates is taken to be 2 arcsec, it 
results in a sample of 219 sources cross-matched. Total fluxes in 
optical and Spitzer bands are c omputed from aperture fluxes as 
described in [Cardamone et al.] ( ,2010j . 
Before selecting sources according to their redshift we must 



choose the bands we will use: the redshift range will be selected 
to ensure a good sampling of the UV range. All the optical broad 
bands of the MUSYC survey are considered (U, U38, B, V, R,I,z) 
and J,H,K bands are added when available. We also add data 
from intermediate band filters whose 5cr depth is fainter than 
25 ABmag (table 2 of Cardamone et al. (2010)). Practically we 
use IA-427,445 ,484,505,527,550,574,598,624,65 1,679 and 738. 
Figure 1 illustrates the rest- frame wavelengths corresponding to 
these filters as a function of z. 

In order to ensure a good photometric sampling around the UV 
bump (2175 A) we select galaxies with redshifts between 0.95 
and 2.2. In this redshift range we have more than 10 photometric 
bands available in the UV (1300-3000 A) and a good sampling 
around 2175 A (Figjlj). 76 galaxies fulfill the redshift condition. 
To ensure that sources at z < 1.6 have data below rest-frame 
1800 A, they must have been detected in the U and U38 filters 
above 5cr (see Table 3 of "Cardamone et al.] ( |201 0)). This leads us 
to 39 galaxies all of which are also strongly detected in the other 
broad bands as well as above 5 cr in the IA-427 filter. When z > 
1.6, measurements beyond the U filter are sufficient, so we select 
sources detected with SNR > 5 in the B band, adding an addi- 
tional constrain of a 5cr detection an intermediate filter (IA-427 
for z < 1.8 and IA-484 for z > 1.8), 13 galaxies are thus added 
to the selection. At the end we have selected 52 galaxies. As a 
further check, we retrieve the HST images at the coordinates of 
these 52 sources using the GOODS cutout service V2.0 provided 
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by the Multimission Archive at STScI (|http ://archive . st sci . edu|) . 
We discard 22 sources with bright neighbours closer than 4 arc- 
sec to the Herschel detection and whose IR measurements are 
uncertain (the FWHM of the PSF at 100 yum is 6.7 arcsec). Half 
of the remaining sample of 30 galaxies has a disk-like morphol- 
ogy, the other half is dominated by compact sources with few 
interacting-peculiar galaxies while 28 have spectroscopic red- 
shifts. The sample of galaxies is presented in Table 1. The HST 
images for each target are available as online materia l . 7 gal axies 
are identified as X-ray emitters by |Cardamone et al. (2008J) they 
are flagged in Table 1 . 

As described earlier in this section, the present study is based on 
galaxies strongly detected from the UV to the far-IR, whose rest- 
frame UV spectrum is sampled by several intermediate-band fil- 
ters. These highly selective criteria lead to a rather small sample 
of objects which we will study in detail in this paper in order 
to derive the parameters of a dust attenuation curve. Our anal- 
ysis is complementary to other studies at intermediate and high 
redshifts which are based on larger samples but they lack far- 
IR data. These studies have concluded that a UV bump must 
be introduced in the attenuation curve but they were not able to 
quantify neither the amplitude of the bump nor the modification 
of the attenuation curve ( .Ilbert eFaTl [20091 [Knek et al.l [20TT] ) 

3. SED fitting with CIGALE 

SED fitting is performed with the CIGALE code (Code 
Investigating GALaxy Emission) developped by |Noll et aL] 
( 2009b ). This is a physically-motivated code which derives prop- 
erties of galaxies by fitting their UV-to-far-IR SEDs. CIGALE 
combines a UV-optical stellar SED and a dust, IR-emitting com- 
ponent and fully conserves the energetic balance between the 
d ust absorbed stellar emiss ion and its e-emission in IR. We refer 
to |Giovannoli et al.] ( |2010| ) for details on the use of CIGALE to 
fit SEDs of actively star forming galaxies observed both in UV 
and far-IR (rest-frame). In the present work, the sampling of the 
UV range is very good so we expect to constrain the dust attenu- 
ation curve. To model the attenuation by dust, the code uses the 
attenuation law of Ca lzetti et al., ( ,2000j as a baseline, and off'ers 
the possibility of varying the steepness of this law and add ing 



NoU et al. 



2009b') for 



a bump centered at 2175 A. We refer to 

a complete description of the dust attenuation prescription. The 

dust attenuation is described as: 



A(^)=^(^'W) + ^.o,r,^.W))h^ 



(1) 



where Ay = 5500 A, k\X) comes from Calzet ti et al 
(Eq.4) and Z)^o,r,^bW) 



( |2000t 



the Lorentzian-like Drude profile com- 
monly used to describe the UV bump ( [Fitzpatrick & Massa 
[T990l|Non etaL||2009a| ) is defined as: 



2^2 



D 



^o,y,E^ 



E^A'r 



(^2 - ^2)2 + ^2^2 



(2) 



The peak amplitude above the continuum, full width half max- 
imum (FWHM) and central wavelength of the bump (E\), y and 
/lo) are free parameters of the codq^ There is no indication for 
variations of the central wavelength of the bump either in the 
Milky Way or in external galaxies ( Noll et aL| |2009a , and ref- 
erences therein), therefore we fix the central wavelength of the 



£b a nd y are related to the ori ginal parameters c^ and 7fm intro- 



bump /lo to 2175 A. The parameter y, indicating the FWHM, can 
vary between 200 and 500 A in order to cover the range of values 
found in previous studies: the FWHM foun d in the Milky Way 



is 437 A ( [Fitzpatrick & Massa 



2007 ) while 



Noll et al. 



(2009a I 



found a FWHM of 274 A for their composite spectrum of galax- 
ies with 1.5 < z < 2.5. After several trade-off's, the amplitude 
of the bump E\^ is chosen between and 2. For comparison the 
corresponding value for the extinction curve of the Milky Way is 
12007) ), 1.63 for the LMC2 supersheff 
between 0.5 and 1 for z ~ 2 galaxies 



3.52 (Fitzpatrick & Massa 
(Gordon et al. , 2003 ), and 



(NoU etal, 2009a). 

Finally we choose d between -0.3 and 0.1, this range of values is 
also chosen after several trade-off's, excluding values never cho- 
sen by the code. 

The stellar population created by CIGALE is the combination 
of an old, exponentially decreasing, component (starting at the 
age of the universe at the redshift of the galaxy) and a young 
burst of constant star formation beginning between 0.3 10^ and 
310^ years ago. The e-folding rate of the old stellar population is 
found to be 2 + 0.4 Gyr. Since the UV emission reaches a station- 
ary state after ~ 10^ years of constant star formation, the general 
shape of the UV SEDs analyzed in this work is predominantly 
modified by the eff'ect of dust attenuation ( Calzetti et aH|1994| ). 
The attenuation correction is applied to each component. To take 
into account the diff'erent distributions of stars of diff'erent ages, 
we adopt a visual attenuation Ay for the old stellar component 
equal to half that used for the young component. The exact value 
of this reduction factor (called fan in the code) is very difficult 
to constrain even if a value l ower than 1 is always found ( Noll et 
aL| |2009b|puat et al.|[2QTT] ). In a first step we ran CIGALE with 



fan as a free parameter with values comprised between and 1, 
we obtained a mean value of 0.5 for fan- To reduce the compu- 
tation time we decided to freeze the value of fan, after checking 
that it did not modify the values obtained for the other output 
parameters of the code. 

Dust luminositi es, log(LiR) between 8 and 1000 yum, are com- 
puted by fitting |Dale & He lou ( 2002 ) templates and are linked 
to the attenuated stellar population models (the stellar luminos- 
ity absorbed by the d ust is re-emitted in IR). The validity of 
Dale & Helou p002 ) templates to measure total IR luminosi- 



ties of sources detected by Herschel is confirmed by the studies 
of |Elbaz et al . (2010, 2011 ). Non-thermal sources, not linked to 
the absorbed stellar population can also contribute to the dust 
emission. CIGALE includes dust-enshrouded AGNs by adding 
two PAH- free templates from Siebenmorgen et al?| ( |2004| ) as de- 
scribed in Buat et al. ( 2011| ). These two templates diff'er by the 
amount of silicate absorption and their spectral distribution be- 
yond 20 microns is sim ilar to the mean template constructed by 
MuUaney et al. ( 2011| ). The spectrum below 20 microns is rep- 



duced by Fitzpatrick & Massa (|l990j): E]j - 03/ yl^ and 7 = 7fm 



xAi 



resentative of type 2 AGNs. 
The AGN contribution to Ljr is found to be lower than 10% ex- 
cept for N-17, the only X-ray source of our sample with broad 
emission lines, and N-25 which has the highest X-ray luminosity 
(log(Lx/Lo) = 44.4) and exhibits a power-law spectral energy 
distribution. For these two sources the AGN fraction is found to 
be larger than 50%. The current version of CIGALE is unable to 
fit accurately these SEDs leading to minimum ;^^ values equal to 
16 and 6.4 for N-17 and N-25 respectively. These two sources 
will be discarded when the output parameters of the code will 
be discussed. The case of the other X-ray sources is discussed 
in section 4.2. Examples of full SED fitting are shown in Fig. [2] 
The figures for the full sample are available as online material. 
All the output parameters are estimated through the Bayesian 
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Fig. 2: Examples of SEDs fitted with CIGALE. The best model (x^ minimization) is plotted together with the observed fluxes. 



analysis. Results for the total IR luminosity, Ljr, the dust atten- 
uation at 1530 A, AjjY, and the peak amplitude of the bump, E\y, 
are listed in Table 1 tog ether with their errors, als o evaluated by 
a Bayesian calculation ( Giovannoli et al.| |2Q10| ). The range of 
log(LiR/Lo) is 11.25-12.4 : all the sources are luminous infrared 
galaxies (log(Li^/ Lq) > 11, LIRG) or ultra luminous infrared 
galaxies (log(LiR/Lo) > 12, ULIRG). The mean dust attenua- 
tion is found to be quite large with < Ay >= 0.92 + 0.34 mag 
and < AjjY >= 3.09 ± 1.1 1 mag. Such a large dust attenuation is 
expected for LIRGs and ULIRGs and is favourable to search for 
the imprint of UV bumps in the spectrum of these galaxies. 

4. Evidence for a bump at 2175 A 

4.1. Amplitude of the bump from SED fitting with CIGALE 

Some examples of SEDs in the UV range are shown in Fig. [3] 
The figures for the whole sample are available as online mate- 
rial. These few examples show the excellent sampling of the UV 
SED achieved thanks to the use of medium filters. A dip at 2175 
A is clearly seen in most of the observed SEDs. The best mod- 
els obtained by the x^ minimization are found to reproduce the 



data very well. They all correspond to a bump amplitude larger 
than 0. The median value of Eiy for the best models is 1.5. In 
Fig. |4] the Bayesian estimate of £'b is plotted against the red- 
shift of the sources. All the sources are found to exhibit a bump. 
The uncertainty given by the Bayesian analysis is quite large but 
a null amplitude is always excluded at the la level. No trend 
is found either with the redshift or dust luminosity. The mean 
value < £'b >= 1.26 ± 0.30 corresponds to 35% of the ampli- 
tude of the bump in the average extinction curve of the MW and 
to 76% of that in the LMC2 supershell. (Note that w e are deal- 
ing with attenuation curve and not extinction ones). Noll et al. 
(2009a) measured the amplitude of the UV bump in 68 spec- 
tra of galaxies with 1.5 < z < 2.5 as part of the GMASS sur- 
vey with a formalism similar to that used in the present work. 
They detected robustly a UV bump for 24 % of their sample 
with Eij = 0.9 ±0.1 while Et, = 0.5 ± 0.3 for the remaining 
76% of their sample. Our result is marginally consistent with the 
mean amplitude they measured for their robust detections of a 
UV bump. 

The FWHM of the bump is found to be 356 ± 20A, intermedi- 
ate between the values found in the MW (and the LMC) and in 
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Fig. 3: Spectral energy distributions for A < 3500 A for the same objects as in Fig|2J the units, symbols and lines are the same as in 
Fig p] except for the scale of the x-axis which is now linear. 



galaxies at z ~ 2 (Nolletal. 2009a I. 



The mean slope of the attenuation curve, S, is found to be -0.13+ 
0.12, no trend is found for ^ as a function of £^ as shown in 
Fig|5] Our result is marginally consistent with the Calzetti et al. 
law which corresponds to ^ = 0. 



4.2. X-ray emitters 

Extinction curves of AGN and quasars are found to not display 
the 2175 A feature and the wavelength dependence of the ex- 
tinction in these objects remains controversial between a gray or 
SMC-like extinction ( |Li1|2007[ and references therein). N17 and 
N25, two of the seve n galaxies in our initial sa mple classified 
as X-ray emitters by Cardamone et al. ( 2Q10| ), are discarded 
because of the bad fitting of their SEDs (see discussion in 
section 3). The SED fitting process gives amplitudes of the 
bump for the remaining 5 sources which are diff'erent from 
but in general lower than for non X-ray galaxies (Fig|4]). We can 
quantify this diff'erence: for a Gaussian probability distribution 



of the parameter E^, Eh/ err e,^ > 2 corresponds to E^ > 



with a probability of 95% (errEi, is the standard deviation of 
the distribution). Only 2 out of 5 X-ray galaxies fulfill this 
condition compared to 19 of the 23 non X-ray galaxies (Fig|6]). 
For the first one, N14 (MUSYCID-30046), Et/errE, = 2.7. This 
galaxy has a high X-ray luminosity, log( Lx/LQ) = 44 . 1 with a 
hardness ratio equal to 0.63 (C ardamone et al.[ |2QQ8| ). The fit 
is of average quality with Xmtn - 4.56 and the presence of a 
bump is not obvious from a visual inspection of the SED. It is 
at least partly due to the fact that the best fit model corresponds 
to a large width of the bump: 500 A. The second one is N21 
(MUSYCID-55660). This galaxy (G02:638) is detected in the 
hard X-rays with a signal to noise ratio of 2 by |Giacconi et al. 
(2002) only using Sextractor and not with a method based on 
wavelet transform; [Alexander et al.| ( |2003| ) did not detect this 
source, so the X-ray detection may be spurious. Its SED shows 
a clear absorption feature at 2175 A with Eh/errE}, =5.1 and the 
SED fitting is good (Xmin = 1.71). Given the small number of 
X-ray sources in our sample as well as the uncertainties inherent 
to the two sources for which a bump is measured above 2cr, we 
cannot draw any firm conclusion about the presence or not of a 
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Fig. 5: 6 versus the amplitude of the bump Eij. X-ray emitters are 
encircled with squares 



UV bump in the attenuation curve of X-ray galaxies. 
The slope of the attenuation curve, 6, is found ^ for two X-ray 
sources and ^ -0.2 for the 3 other ones. Therefore a steep 
attenuation curve is favoured, not as steep as the SMC extinction 
curve but excluding a gray extinction which corresponds to 
positive values of d. 



4.3. Average dust attenuation curve 

The resulting mean dust attenuation curve obtained for the whole 
sample is plotted in Fig. |7] All parameters are calculated by a 
Bayesian analysis. The average dust attenuation curve we obtain 
can be written: 



Ay ~ 4.05 \Av 



-0.13 



where Ay = 5500 A, k\A) is given in |Calzetti etaLl ( |2000| (Eq.4) 
and 



D 



^0,7, Eh 



W) 



1.26x3562^2 



(^2-21752)2 + ^2x3562 



(4) 



(3) 



Strictly speaking this average attenuation curve applies to 
the young stellar component only (as explained in section 3), 
the dust attenuation applied to the old component being reduced 
by a factor 2. In practice, the SEDs of our galaxies appear to 
be dominated by the young stellar component (modeled as a 
constant star formation rate over ~ 10^ years), especially in UV. 
As a result the global attenuation curve is found to be quite 
similar to that of the young stellar component alone. 

The uncertainty on the average value of 6 implies quite 
a large dispersion of the attenuation curve at A < 2500 A 
(Fig.[7| and of the determination of dust attenuation at UV wave- 
lengths. For a visual attenuation Ay = I mag the difference be- 
tween the dust attenuation at 1500 A obtained by using Eq. 3 
(S = -0.13) or a Calzetti et al. law (S = 0) reaches 0.5 mag. A 
value of ^ 9^ also implies a change in the value of the e ffec- 
tive total obscuration Ry = Ay /Eb-v. [Calzetti et al.| ( |200Ql ) ob- 
tain Ry = 4.05 + 0.80 by comparing the predicted and observed 
absorbed emission (i.e. by performing an energetic budget). 
Adopting a value of ^ = -0.13 leads to Ry = Ay/E e-y = 3.7, 
still consistent with the value found by [Calzetti et al.| ([2000). 
The introduction of a bump of moderate amplitude in the dust 
attenuation curve (Eq. 3) leads to an increase of A(A)/Ay at 
^ 2000 A of ^ 0.35 mag as compared to the same curve with 

DAo,y,EM = 0. 

We can compare our average du st attenuation cu rve to other 
attenuation curves in the literature. 'Chariot & Fall (2000]) pro- 
pose a recipe consisting in attenuating the stellar emission by a 
factor proportional to A~^-^ , and in reducing the normalization of 
the attenuation for stars older than 10^ years; the resulting effec- 
tive absorption curve they obtain for a constant star formation 
rate over 3 10^ years (their Fig. 5) is reported in Fig.JTJ A period 
of 3 10^ years is close to the mean duration of the burst of star 
formation we obtain for the young stellar population (10^ years, 
cf. section 3). Apart from the presence of the bump, not consid- 
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ered b y [Chariot & Fall|(|2000 ), our mean attenuation curve and 
that of Chariot & Fall '{ 200()] rexhibit similar shapes and steepen- 
ing for decreasing v^avelengths. 

Although we are dealing with attenuation and not extinction, 
we also compare our resulting curve with the extinction curves 
found in the Milky Way and the Magellanic Clouds. The rela- 
tive extinctions are compared in Fig. [8] The general shape of the 
extinction curve of the LMC2 supershell is consistent with our 
attenuation curve, whereas the MW and SMC extinction curve 
show flatter and steeper variation, respectively, at short wave- 
lengths. 

The exact shape of the extinction curve (slope in the UV range 
and amplitude of the bump) cannot be constrained with the 
present data since the attenuation curve we derive strongly de- 
pends on the distribution of both stars and dust. Witt & Gordon 
( |2000[ ) show that the shape of the attenuation function also de- 
pends on the total amount of dust attenuation. We do not observe 
any variation of E\y and S with Ayv- it can be explained by the 
homogeneity of our sample in terms of dust attenuation associ- 
ated to the large uncertainties on the determination of these pa- 
rameters. Radiative transfer models in opaque disks with various 
configurations of stars and dust (especially in a clumpy medium) 
are able to reduce the amplitude of the UV bump in the eff'ective 
attenuation with a dust composition similar to that of the MW or 
the average LMC, the attenuation curve becoming grayer when 
the opacity increases (e.g. [Witt & Gordon |2QQQi Pierini et al .| 
[20041 [Inoue et al.[ [20061 [Panuzzo et al.l [2Q07| ). These models 
also predict a flattening of the attenuation curve in UV (as an ef- 
fect of the gray attenuation) as compared to the original extinc- 
tion curve. The average attenuation curve we deduce increases 
steeply at short wavelengths: it seems to favour the scenario of 
a deficiency of bump carriers with respect to the case of MW or 
average LMC. We will see in sectio n 5.2 that our source s roughly 
follow the Ajjy-p relation found by [Meurer et al.[ ( p^99] ) for local 
starbursts which is wel l explained by a clumpy shell geometry 
dWitt & Gordon[ ( |20Q0| ); [Calzetti[ ( [2QOT| ). Such a dust/star con- 
figuration favours the interpretation of a moderate UV bump in 
the extinction curve itself in agreement with our conclusion (e.g. 
[Witt & Gord on, 2000; NoU et al. , 2009a). 

5. Slope of the UV continuum from broad band 
photometry and dust attenuation 

5. 1 . Determination of the UV slope 

5.1 .1 . Influence of the bump 

With the imprint of the UV bump clearly visible in our sam- 
ple, the assumption of a simple power-law for the entire UV 
continuum is no longer valid, and its influence in the mea- 
sured slopes can be quantified. A single power law might be 
valid if we exclude the spectral region of the bump (Fig [3]). 
Since our fluxes are defined in Jansky we define Of as /y oc A^ 
with a = p -\- 2 where fi corresponds to the original definition 
/i(erg cm"^s"^A"^) oc A^. From Fig. [sjit is clear that we can 
extend the model of a power-law up to 3000 A. Therefore we 
calculate the slope of the UV continuum in several ways. First 
we use all the data between 1200 and 3000 A, avoiding the bump 
range (1975-2375 A) and perform a linear regression. The result- 
ing value of a will be taken as our reference a^ef. Then we define 
3 wavelength ranges 1500-2000 A (Bl) , 1500-2500 A (B2) and 
2500-3000 A (B3) which mimic broad band filters and combine 
two of these bands to estimate a. The definition of these 3 bands 
is illustrated in Fig|9]for a galaxy exhibiting a strong bump. B2 
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Fig. 7: Average dust attenuation curve. The mean curve found in 
this work is plotted as a thick line. The thin lines correspond to 
variations of ± Icr of the amplitude of the b ump, E^ and the slope 
of the continuum 6; see text for details. The |Calzetti et al.| ( [2000| ) 
attenuatio n curve (dotted blue line) and the effective absorption 
curve of Char iot & Fall ( 2000*) for a starburst age of 3 10^ years 
(dot-dashed red line) are shown on the figure for comparison. 




Fig. 8: E(A - V)/E(B - V). The mean curve found in this work 
is plotted as a thick line together with the extinction curve of 
the MW (red dotted line) from |Whittet| ( [2003] ), LMC2 supersheU 
(green short dashed line) and SMC (blue dot-dashed line) from 
[Gordon etal.| ( |2003| ). 
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Fig. 9: Definition of 3 bands for ID 2733. The best model found 
for this galaxy is plotted and the range of wavelengths corre- 
sponding to B 1, B2 and B3 is in blue, green and red respectively. 
The black dots are the observed fluxes. 



overlaps the bump, therefore we expect some efl'ects when the 
slope is calculated by considering it. The slopes of th e UV con- 
tinuum are calculated following the prescription of jKong et al. 
( |2Q04| ) adapted to our notation: 



a = 



logifvjfvj) 



(5) 



where /y. and /y . are the mean flux densities per unit frequency 
through the filters i and j whose eff'ective wavelengths are At and 
Aj. In FigfTOlthe a values obtained by combining Bl and B2 
(0^12), B2 ana B3 (0^23) and Bl and B3 (^13) are compared to 
Qfref measured on the entire UV SED. As expected ai3, calcu- 
lated avoiding the bump is found to correlate very well with a^Qf 
(Fig 10 top panel, R = 0.92). The combination of B1-B2 leads 
to bluer slopes whereas the slopes calculated with B2 and B3 
are redder than a^Qf. These eff'ects are expected since the pres- 
ence of the bump increases the dust attenuation in the B2 band. 
The amplitude of the shifts are < au - c^ref >= -0.46 ± 0.26 and 
< a23-arQf >= 0.28+0.40 whereas < ^is-^ref >= -0.09+0.25. 

5.1.2. GALEX-like UV slope 

The FUV and NUV bands of GALEX are commonly used to 
mea sure the slope of the UV continuum in the nearby universe 



(e.g. Seibert et al. 



(Morrissey et al.| 



2005]). Since the NUV band (1771-2831 A) 
2005 1 ) overlaps the UV bump we might ex- 



pect some underestimation of the slopes, ckgalex is an output of 
CIGALE: the spectral energy distribution of each model gen- 
erated by the code is integrated in the FUV and NUV (rest- 
frame) filter bands and ^galex calculated with the above for- 
mula. Then the expected value of ^qalex for each g alaxy of our 
sample is estimated through a Bayesian calculation ("Giov annoK 
et al.[ [2010). The correlation between q^galex and a^Qf is very 
good (R = 0.91). As expected ckgalex is typically lower (bluer) 
than Qfref but the systematic diff'erence remains small with < 
<^GALEX - o^ref >= -0.19 + 0.24. This moderate shift, smaller than 
the one found considering 0^12, is probably d ue to the wi de large 
bandpass of the NUV filter (1060 A, ( Morrissey et al.| [2005)), 
much wider than the bump itself. Therefore, as long as galaxies 
exhibit bumps in their attenuation curve similar to those found 





Fig. 10: Various estimates of a compared to the reference value 
Qfref measured on the full SED , top: 0^13 versus a^Qf', middle: 
au (blue empty circles) and 0^23 (red empty squares) versus Ofref ; 
bottom: o^galex calculated with the rest-frame FUV and NUV 
fluxes versus a^ref • On each figure, the dashed lines represent the 
1:1 line, the solid lines are the results of linear regressions. 



in this study (1/3 that of the Milky Way), the slopes measured 
with the GALEX data are found reliable. Burgarell a et al. (2005j ) 
showed that Milky Way-like bumps have a large influence on the 
determination of the UV slope with GALEX filters. 
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Fig. 1 1 : UV dust attenuations from CIGALE plotted against the 
slope of the UV continuum calculated with the full SED. The to- 
tal IR luminosity is colour-coded. The solid and dashed lines cor- 
respond to the inner and total relations of |Overzier et al.| pQll| ) 
respectively. The error bars come from CIGALE (Bayesian esti- 
mates of the parameters and their associated errors). The red line 
is the result of the linear regression 



5.1 .3. Impact of photometric errors 

It is worth noting that the determination of the UV slope through 
two broad band filters is very sensitive to photometric errors as 
well as to the interval between the eff'ective wavelengths of both 
filters i.e. the factor log(/l///ly) in Eq.(3). We can re-write Eq.(3) 
in terms of magnitudes instead of fluxes: 



a=/5 + 2 = -0.4 



(nii - ruj) 



(6) 



Using the FUV and NUV filters of GALEX the relation becomes 
a = -2.22 (rrii - rrij) and a typical error of 0.1 mag in FUV and 
NUV translates to an error of 0.3 on a. It is obvious from Eq.6 
that increasing the wavelength diff'erence between the two fil- 
ters reduces the uncertainty on the determination of a and/?. For 
the domain of validity of the power law adopted in this work 
(A < 3000A, cf. section 5.1.1), using filters with the largest pos- 
sible interval between rest-frame wavelength, i.e. 1300 A (to 
avoid the flattening of the attenuation cu rve for A < 1200 A 



(Leitherer et al. 2002 Buat et al. 20021 and 3000 A reduces 



the error by a factor ~ 2 (this time or = -1.1 l(m/ - rrij)) as com- 
pared to the use of GALEX filters. 

To summarize, given the uncertainty on the presence and 
strength of a UV bump in a galaxy, one must avoid as much 
as possible the region aff'ected by the bump when measuring the 
UV slope using broad band filters and use filters with eff'ective 
wavelengths that are as widely separated as possible from each 
other. If the region of the bump cannot be avoided, the attenua- 
tion curve presented in Eq.3 might be used to interpret the data. 

5.2. Dust attenuation corrections 



[Meurer et al.| ( |1999| ) first deduced a relation between Ajjy and 
13 for local starbursts, avoiding the bump area on the spectra ob- 
tained with lUE. Sev eral calibrations have been proposed since 
this pioneering work ( [Calzetti et aL| |2000t [Seibert et aL) |2005| 




1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 
^UV,CIGALE ("^^g) 



Fig. 12: Comparison of UV dust attenuations: the values of Ayy 
directly obtained from the SED fitting with CIGALE are plotted 
on the X-axis and the estimates of Auv from diff'erent UV slopes 
and Eq. 7 on the y-axis . 



Qverzier et al.| 201 1| ). One issue of the preliminary analyses was 
that the UV data came from lUE whose aperture sampled only 
the inner part of the galaxies. Qverzier et a l. ( 201 1) re-investigate 
the topic by using GALEX images defining an 'inner' and 'to- 
tal' relation, the first one based on UV fluxes measured into an 
lUE-like aperture , the second one corresponding to the total UV 
fluxes. These relations are: 

Auv = 4.54 -h 2.07/5 = 4.54 -h 2.07(a^ - 2) 

for the inner relation and 

Auv = 3.85 -h 1.96/5 = 3.85 -h 1.96(Qf - 2) 

for the total relation. In Fig[TT] these relations are compared 
to the values of a^Qf and Auv,cigale (the dust attenuation ob- 
tained from the SED fitting process and a Bayesian calculation). 
Whereas LIRGs seem to follow more or less the local starburst 
relations, the ULIRGs are found above the local relations which 
is in agreement with p r evious studies at low and high redshift 
(Takeuchi et al.l |20T0l |HoweU et al.| [2QTq1 [BurgareUa et al. 



2007| |Reddy et al.| |2008| ). |Buat et al.| ( |201QJ ) performed a se- 
lection at 250 fim at z < 0.3 and found that their sample (in- 
cluding LIRGs) lie below the local starburst relation. This dif- 
ferent behavior is likely to be due to a selection eff'ect biased 
towards more quiescent galaxies when selected at 250yum than 
in the mid-IR or around 60 jim and to the use of the GALEX 
colour to measure the UV slope at z > 0. 
A linear regression on the whole sample gives: 



Auv,ciGALE = 1.46 (±0.21) aref + 1.50 (±0.29) 



(7) 



with a correlation coeflficient R = 0.79. The correlation is rather 
weak resulting in a large dispersion with a RMS scatter of 0.73 
mag. This error directly aff'ects any estimate of Ljr based on re- 
lations between Ljr/Luv, a quantity tightly related to Auv (Buatj 



et al. 2011 and references therein), and the slope of the UV 
continuum. 

We can also compare dust attenuation corrections obtained 
by applying Eq. 7 to the various estimates of a defined in sec- 
tion 5.1.1 and 5.1.2. The results are reported in Fig. 12 and com- 
pared to the dust attenuation given by CIGALE: Auv,cigale. A 
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large dispersion is found, as a consequence of the loose rela- 
tion found initially between Ajjy and Ofref (Fig H ): < AuY,aM ~ 
^uv,ciGALE >= -0.005 ± 0.707 mag. When using 0^13 the sys- 
tematic shift remains moderate < Auv,ai3 ~ ^f/v,ciGALE >= 
-0.18 ± 0.58 mag. ckgalex, ^12 and CK23 give worse estimates 
with substancial shifts and large dispersions: < Auv,aGALEx ~ 

A[/y,ciGALE >= "0.32 ± 0.58 mag, < Auy^^u - ^UV,CIGALE > = 

-0.73 ± 0.70 mag, < Auv,c.23 " ^uv,cigale >= 0.37 ± 0.71 mag. 
Therefore, using UV slopes derived from broad band colours 
leads to large uncertainties in the derivation of dust attenuation 
(~ 1 mag) even if the bump area is avoided and good photometry 
is available. This is due to the intrinsic dispersion found between 
the slope of the UV continuum and the amount of dust attenu- 
ation. The presence of a bump in one of the filters may add a 
systematic shift which can reach 0.7 mag for the galaxies of our 
sample which is likely to be representative of star forming galax- 
ies at redshifts between 1 and 2. These results hold for galaxies 
with a moderate bump in their dust attenuation curve (35% of the 
mean amplitude found in the MW) and we expect larger eff'ects 
for larger bumps. 



6. Conclusions 

We have identified the presence of a UV bump in a sample of 
galaxies of the CDFS at a redshifts between 0.9 < z < 2.2, 
selected to be observed in intermediate and broad band optical 
filters as well as at mid-and far-IR wavelengths by Spitzer 
and Herschel. This study demonstrates the capabilities of 
intermediate-band photometry to give details on SED of galax- 
ies, and to derive physical properties. 

The average dust attenuation curve we deduce is well described 
by a modified Calzetti et al. law slightly steeper than the original 
one and with a UV bump at 2175 A whose amplitude is ~ 
35% of the MW one. We propose an analytical expression of 
the average attenuation curve which can be used to correct the 
SEDs of galaxies for dust attenuation. The moderate amplitude 
of the bump together with the substantial slope of our average 
attenuation curve in the UV argue for a deficit of UV bump 
carriers in our sample galaxies as compared to the MW. 
Our sample contains seven X-ray galaxies. Five of them are 
reasonably fitted: a steep attenuation curve is favoured excluding 
a gray extinction. The sample is too small to draw any firm 
conclusion about the presence of a UV bump. In any case the 
amplitude of the bump, if any, is smaller than that found for non 
X-ray galaxies. 

The presence of a bump in the mean dust attenuation curve 
has implications for the derivation of the slope of the UV 
(rest- frame) continuum of galaxies from broad band data alone. 
The power law modeling of the UV spectral distribution of our 
galaxies can be extended up to 3000 A but the bump area must 
be avoided. When a broad band filter overlaps the bump, despite 
its moderate amplitude, the error on the determination of the 
slope may reach -0.5. The use of GALEX filters leads to an 
underestimation of the slope of only 0.2, even if the NUV filter 
overlaps the bump. This is likely to be due to the very large 
bandpass of this filter. However, the proximity of the two central 
wavelengths of the GALEX filters implies a large uncertainty 
in the determination of the slope. It is stressed that a larger 
wavelength baseline would give a better measure. 
Dust attenuation estimated with UV slopes, even when the latter 
are reliable, remains uncertain, with an RMS error of 0.7 mag. 
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Fig. 1: Spectral energy distributions for the full sample and HST 
images obtained from the GOODS cutout service V0.2. The best 
fit is represented with a solid line. The rest- frame wavelength is 
plotted on the x-axis and the fluxes in Jy on the y-axis. 



Fig. 1: cont'd 



V. Buat et al.: GOODS-Herschel: evidence for a UV extinction bump in galaxies at z > 1, Online Material p 2 




N=10 z=1.04 




2000 3000 

A(A) 




B ■ - ■ ■ 




I ■ ■ ■ 


z ' ' ■ 







10 












i 

1 


B 




V 




I 


z 









.f» 




-f> 


^ 




m 


-7 

















N = 8 z=1.02 



2000 3000 

A(A) 




N=ll z=1.09 




2000 3000 

A(A) 




11 



B 

• 
f 


V 

1 


I 

• 

1 





B ■ ' ' 


V ' ■ ' 




z ■ 




N=12 z=1.10 




2000 3000 

A(A) 



3 4 5 6 7 

A(A) 



12 







I z 



B ■ ' ' V ■ ■ ■ 


I ■ ■ ■ 


z ■ ■ ■ 



Fig. 1: cont'd 



Fig. 1: cont'd 



V. Buat et al.: GOODS-Herschel: evidence for a UV extinction bump in galaxies at z > 1, Online Material p 3 




N=16 z=1.10 




2 - 



2000 3000 

A(A) 



5 
A(A) 



13 



16 



B ' ' ■ 

« 


V 


I " 

• 


z ■ ■ ■ 



B ■ ■ ■ V ■ ■ ■ 
# 


I ■ ■ ■ 


z ■ ■ ■ 

• 




2000 3000 

A(A) 




2000 3000 3 

A(A) 







14 














2 
1 


B ■ 




V 




I 




z 




1 

-1 








*♦- 




^ 




^ 


-7 



















17 



B V I Z 




2000 3000 

A (A) 



N=18 z=1.22 



2000 3000 

A(A) 




15 



B ' 




V ' 




I 




z 






- 




m 




^ 




» 




^ 




, 











1 
-1 

-2 



B ■ ■ ■ :V ■ ■ ■ 


I ■ ■ ■ 


z ' ■ ■ 



Fig. 1: cont'd 



Fig. 1: cont'd 



V. Buat et al.: GOODS-Herschel: evidence for a UV extinction bump in galaxies at z > 1, Online Material p 4 




N = 22 z=1.38^ 




2000 3000 

A(A) 




19 



B 

4 


V ■ ■ ■ 

* 


I ■ ' ' 

* 


z ' ' ' 


m 







22 














I 

1 


B 




V 




I 




z 









r 




t 




# 




# 


-1 

-7 




















N = 23 z=1.38 




2000 3000 

A(A) 








20 








2 


B 




V 


I 


z 


1 




«te' 


m^ 


411^ ^ 


m - 


-1 












-7 













23 



B ' ' ' 


V ' ' ' 

v^ - 


I ■ ■ ■ 


z ■ ■ ■ 




-2.5 - N-24 z-1.52 




2000 3000 

A(A) 



2000 3000 

A(A) 




21 



24 



B ■ '^ ■ 


• 


I t ' , 


z ■* ■ 





V 


I ■ ' ■ 


z ' • ■ 



Fig. 1: cont'd 



Fig. 1: cont'd 



V. Buat et al.: GOODS-Herschel: evidence for a UV extinction bump in galaxies at z > 1, Online Material p 5 




N = 28 z=1.99 




2000 3000 

A(A) 








25 














2 


.B 




-v 




,1 




,z 




1 























• 




•- 




« 




• 


-1 


















-7 






. 








, 





28 



B ■ ■ ■ 


V 




I ■ ■ ■ 
1^ 


z 




N-29 z-2.00 




2000 3000 

A(A) 




26 



29 



B " 


V ' ' ' 


f 


Z 

»- 



B 


V 


I ■ 


z ' ' ' 

M 




2000 3000 

A(A) 



N = 30 z = 2.13 




2000 3000 

A(A) 




27 



B ■ 


> 


V 






,«-i55J 


.^^^,^!-A.^;^.^A:. W,v! 



30 



:B " ■ ■ V ■ ■ ■ 

* * 
* • 


I ■ ' ■ 


z ■ ■ ■ 1 

•- 
• 



Fig. 1: cont'd 



Fig. 1: cont'd 



